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ABSTRACT 


Tornado  Identification  from  Analyses  of  Digital  Radar  Data 

(Oecetnber  1976) 

Donald  Wayne  Pittman,  B.S.,  St  Louis  University 
Directed  by:  Dr.  George  L.  Huebner 

An  Investigation  was  conducted  to  determine  whether  tornadoes 
presented  a unique  signature  In  analyses  of  digital  radar  data  from 
central  Oklahoma  during  the  Spring.  The  data  were  collected  by  the 
ID-cm  WSR-57  radar  at  the  National  Severe  Storms  Laboratory  at  Norman, 
Oklahoma.  Three  types  of  numerical  analyses  were  used  In  this  study: 
constant-altitude  reflectivity  maps  (CAZN),  total  vertically-summed 
reflectivity  maps  (TVSZ),  and  partial  vertically-summed  reflectivity 
maps  (PVSZ),  with  greatest  emphasis  placed  on  the  PVSZ  maps.  Presen- 
tations covering  a 100-km  square  were  constructed  at  either  5 or  10 
min  Intervals. 

From  the  analysis  of  three  case  studies  during  1974  and  1975, 
that  contained  five  tornadoes.  It  was  concluded  that  tornadoes  did  not 
produce  a singular  Identifying  signature  In  analyses  of  digital  radar 
data,  but  rather  produced  a combination  of  features  which  Indicated, 
with  a high  probability,  the  presence  of  a tornado.  Such  features 
were  the  appearance  of  a small  area  of  reduced  reflectivity  known  as  a 
Bounded  Weak-Echo  Region  (BWER),  a tilt  of  the  core  of  the  storm  toward 
the  BU. ' a rapid  decrease  In  the  upper-level  mass  of  the  storm  as 
1nd1ca». > ^ rapid  decrease  In  the  reflectivity  of  the  upper  PVSZ. 
Further  Investigation  Is  required  before  any  definitive  value  of  prob- 
ability can  be  assigned  to  this  Identification  technique. 
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CHAPTER  I 
INTRODUCTION 

Severe  local  storms  are  composed  of  Intense  thunderstorms  that 
produce  heavy  rain,  structurally-damaging  winds,  crop-damaging  hall, 
tornadoes,  or  any  combination  of  the  latter  three.  Each  year  through- 
out the  world,  these  types  of  storms  cause  millions  of  dollars  of  dam- 
age and  the  loss  of  numerous  lives.  In  1975  In  the  United  States  alone, 
severe  storms  were  Involved  In  the  death  of  60  persons,  and  property 
damage  for  a single  storm  reached  a new  high  of  160  million  dollars. 

It  then  Is  no  wonder  that  researchers  are  continually  at  work  to  try 
to  find  better  methods  of  detecting  and  forecasting  the  occurrence  of 
these  storm. 

The  Need  for  this  Investigation 

While  better  methods  of  detecting  and  forecasting  severe  storms 
would  produce  more  accurate  warnings,  their  benefit  to  the  populace 
would  be  only  as  great  as  the  response  generated.  To  date,  the  accu- 
racy of  severe  storm  warnings  has  been  such  that  populace  response  has 
remained  minimal  [Muench , 1976].  Commercial  businesses  must  weigh  the 
cost  of  preparation  and  action  during  each  false  alarm  against  the  risk 
of  destruction.  Hence,  due  to  the  relatively  Infrequent  occurrence  of 
severe  storms  at  a given  location,  the  cost  of  actions  taken,  and  the 
Inaccuracy  of  warnings,  the  response  of  businesses  also  has  remained 


The  citations  on  the  following  pages  follow  the  style  of  the  Jour- 
nal of  Geophysical  Research. 
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minimal.  Therefore,  In  order  to  Improve  the  response  and  minimize  the 
loss  of  life  and  property,  the  accuracy  of  severe  local  storm  warnings 
must  be  Improved.  It  Is  with  this  fact  In  mind  that  this  Investigation 
Is  conducted. 

The  Scope  and  Objectives  of  this  Investigation 

This  Investigation  examines  a new  method  of  detecting  the  most 
damaging  severe-storm  phenomenon,  the  tornado.  Specifically,  the  ob> 
jective  of  this  study  has  been  to  determine  If  storms  producing  a tor- 
nado or  funnel  cloud  present  a unique  signature  In  analyses  of  digital 
radar  data.  Information  recently  presented  seems  to  Indicate  that 
there  Is  a high  correlation  between  the  occurrence  of  tornadoes  In  Ok- 
lahoma In  the  Spring  and  the  appearance  of  a small  transient  area  of 
reduced  reflectivity  known  as  a bounded  weak-echo  region  (BWER)  and  an 
acute  tilt  of  the  core  of  maximum  reflectivity  [Canipe,  1973].  If 
this  Information  Is  correct,  the  meteorologist  would  be  provided  with 
a numerical  tool  for  Identifying  small  areas  In  which  a tornado  Is  oc-  l 

I 

curring. 


To  evaluate  this  hypothesis,  multi-tilt  digital  radar  data  were  I 

manipulated  by  computer  to  produce  several  types  of  numerical  displays.  | 

These  displays  were  developed  for  times  during  which  severe  weather  was  ' 

forecast  to  occur,  not  just  during  the  times  a tornado  was  observed  to  ; 

occur.  In  this  manner  an  objective  evaluation  of  the  correlative  . | 

points  could  be  conducted.  It  Is  not  the  Intent  of  this  Investigation  j 

to  complete  a statistical  analysis  of  the  accuracy  of  this  observa- 
tional tool,  but  rather  to  contribute  case  studies  toward  this  goal. 


Present  Status  of  the  Problem 


In  general,  the  task  of  studying  severe  storms  Is  a difficult  one 
and  progress  has  been  slow.  The  complexity  of  severe  storms  Is  such 
that  modeling  techniques  have  been  Ineffective.  Hence,  the  Investiga- 
tor must  rely  upon  nature  to  produce  the  dreaded  phenomenon  In  a timely 
manner  at  an  observable  location  In  order  to  gain  Insight  Into  Its 
structure  and  dynamics. 

Invariably,  the  study  of  severe  storms  Involves  the  use  of  mea- 
suring equipment;  however,  the  equipment  tends  to  measure  the  total 
effect  of  many  simultaneously-occurring  physical  processes,  thereby 
leaving  the  Investigator  at  a loss  as  to  which  Is  causing  what. 

This  does  not  mean  that  all  Instrunentatlon  Is  of  little  value, 
for  significant  progress  has  been  made  In  certain  areas  of  severe-storm 
Identification  and  forecasting,  even  though  the  Internal  processes  have 
not  been  thoroughly  understood.  Some  of  these  areas  will  be  examined 
next. 

The  Weather  Radar 

Over  the  last  20  yr,  the  most  widely-used  tool  for  severe-storm 
Identification  has  been  radar.  Even  though  early  studies  concentrated 
on  the  morphology  and  travel  of  radar  echoes,  researcher*  did  determine 
that  storms  producing  severe  weather  exhibited  certain  unique  charac- 
teristics. Among  the  more  Important  were  the  following:  that  storms  of 
greater  severity  often  moved  with  a component  to  the  right  of  the  path 
of  smaller,  less  Intense  storms  [Battan.  1959];  that  storms  of  greater 
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severity  had  higher  echo  tops  [Battan,  1959];  that  tornadoes  were  asso- 
ciated with  echoes  which  were  qualitatively  very  Intense  and  which 
sometimes  displayed  a hook  or  "6''-shaped  appendage  [Bigler,  1955];  and 
that  tornadoes  tended  to  develop  on  the  right  rear  side  of  the  Intense 
parent  echo  and  not  beneath  It  [Penn  et  al.,  1955]. 

The  difficulty  with  these  characteristics  was  that  they  were  quail* 
tative  rather  than  quantitative.  Further,  although  the  majority  of 
severe  storms  tended  to  follow  this  set  of  characteristics,  a sizeable 
number  did  not.  Donaldson  [1965]  attributed  part  of  this  latter  diffi- 
culty to  the  geographical  variation  of  the  characteristics.  He  found 
that  storms  in  New  England  that  produced  tornadoes  and/or  large  hall 
were  characterized  by  echo  heights  which  were  significantly  greater 
than  the  average  for  that  area,  yet  no  such  anomaly  was  onserved  In 
severe  storms  of  the  Southwest.  Thus,  Investigators  began  restricting 
their  efforts  to  certain  geographical  locations.  Although  the  number 
of  wel 1 -documented  storms  In  a geographical  area  was  relatively  small, 
researchers  did  have  sufficient  Information  to  determine  that  the  num- 
ber of  exceptions  to  any  set  of  Identifying  characteristics  was  still 
significantly  greater  than  zero. 

As  no  single  total ly-rellable  Identifier  of  severe  weather  had 
been  found,  and  It  did  not  appear  that  one  was  going  to  be  found,  re- 
searchers began  to  conduct  statistical  analyses  of  various  radar  param- 
eters to  determine  which  set  of  characteristics  at  a given  location  was 
most  reliable  In  Identifying  severe  local  storms.  The  results  were 
surprising.  Using  only  the  maximum,  equivalent,  radar-reflectivity 
factor  (Z^),  Boyd  and  Husll  [1970]  found  that  they  could  Identify  90 
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per  cent  of  all  hailstorms  In  western  Nebraska  while  keeping  their 
false-alarm  rate  to  only  15  per  cent.  Wilson  [1971]  found  that  the 
best  Indicator  for  severe  weather  In  Oklahoma  was  the  size  of  the  con> 

5 

tiguous  area  of  severe  weather  having  a radar  reflectivity  of  3 x 10 
mm^m"^.  Brandes  [1972]  found  that  Boyd  and  Musll *s  technique  also 
worked  In  Oklahoma,  but  that  the  critical  value  for  the  maximum  was 
greater  than  for  Nebraska.  One  common  difficulty  among  these  studies 
was  that  each  used  only  the  0-deg  elevation  scan  angle,  generally  be- 
cause of  the  difficulty  In  collecting  and  processing  multi-tilt  data. 

Digital  Radar  Data 

To  overcome  the  problem  of  collection,  storage,  and  processing  of 
multi-tilt  data,  researchers  turned  to  digitizers,  magnetic  recorders, 
and  computers.  The  digitizing  units,  called  Digital  Video  Integrator 
Processors,  took  the  returned  signal  from  radar  targets,  quantified  It 
Into  discrete  levels,  and  computed  an  Integrated  Intensity  level  for  a 
specific  area.  The  result  was  then  stored  on  magnetic  tape  by  a digi- 
tal recorder,  the  data  being  retrievable  by  means  of  a spherical  re- 
ference system  of  range,  azimuth  angle,  and  antenna  tilt. 

Digital  radar  data  readily  lend  themselves  to  computer  manipula- 
tion. Greene  [1971]  developed  a numerical  method  to  display  values  of 
Z-  on  constant-altitude  reflectivity  maps  (CAZM)  similar  to  the  con- 
stant-altitude  plan-positlon-indicator  (CAPPI)  presentations  of  Mar- 
shall [1957].  Greene's  method  converted  the  stored  spherical  data  to 
a three-dimensional  rectangular  presentation  by  the  use  of  a quadratic 
Interpolation  scheme.  Such  a presentation  Is  desirable  for  vertical 


I 

1 

1 


i 

i 


5 


\t 


!*■ 


6 


and  horizontal  Integration  and  sumnatlon  techniques. 

Greene  developed  one  such  analytical  tool  by  converting  reflec- 
tivities to  liquid-water  content,  M,  and  then  vertically  Integrating  M 
through  the  depth  of  the  storm,  thus  yielding  the  vertically-integrated 
llquld-^ater  (VIL)  content.  He  found  that  tornado -producing  systems 
were  characterized  by  high  values  of  VIL  for  approximately  1 hr  prior 
to  the  tornado  and  by  a rapid  Increase  to  even  greater  values  during 
the  development  of  the  confirmed  tornado.  Later,  however,  other  re- 
searchers such  as  Vogel  [1973]  and  Canipe  [1973]  found  that  storms 
other  than  tornado-producers  also  were  characterized  by  high  values  of 
VIL.  In  his  Investigations,  Canipe  [1972,1973]  examined  other  possible 
uses  of  the  CAZM,  Including  the  horizontal  Integration  of  liquid-water 
(HIL)  and  the  partial  vertical  Integration  of  liquid-water  (PVIL). 

The  Basis  for  the  Investigation 

During  his  latter  Investigation,  Canipe  found  nine  tornado-pro- 
ducing systems  that  were  characterized  by  the  presence  of  a small 
transient  area  of  reduced  reflectivity,  the  BWER.  It  appeared  In  the 
parent  cell  or  a nearby  satellite  several  minutes  before  tornado  for- 
mation and  dissipated  shortly  before  the  tornado  visibly  disappeared. 
The  region  generally  appeared  first  at  mid-levels  of  the  storm  and 
developed  downward  with  time.  Perhaps  even  more  striking  was  the  fact 
that  no  such  region  developed  In  ten  other  cases  which  produced  heavy 
rain  and  hall  but  not  tornadoes. 

The  weak-echo  region  Is  not  new.  Its  presence  having  been  found 
by  Browning  and  Ludlam  [1962]  during  their  comprehensive  study  of  a 


storm  near  Workingham,  England,  In  July  1959;  however,  virtually  all 
research  Into  this  phenomenon  has  been  conducted  In  an  attempt  to  under- 
stand the  dynamics  of  thunderstorms  rather  than  to  Identify  severe 
storms.  Analyses  of  storms  by  Browning  and  Ludlam,  Marwitz  and  Berry 
[1971],  and  Rinehart  et  al . [1975]  In  areas  of  England,  western  Canada, 
and  the  Upper  High  Plains  of  the  United  States,  respectively,  reveal 
BWERs  that  are  not  associated  with  tornadoes.  Hence,  the  area  of  appli- 
cability of  Canipe's  technique  must  be  limited  to  Oklahoma  until  further 
research  Is  conducted  In  other  southern  locations.  Since  Canj^e  com- 
pleted his  work,  further  confirmation  of  the  association  of  BWERs  with 
tornadoes  In  Oklahoma  has  come  from  studies  by  Bensch  [1975],  Lemon  et 
al..  [1975],  and  McCarthy  et  al . [1975].  Also,  an  analysis  of  a toma- 
do-less  hailstorm  by  Brandes  [1975]  revealed  no  BMER. 

The  terms  "Echo-Free  Vault,"  "Weak-Echo  Region,"  and  "Bounded 
Ueak-Echo  Region"  are  used  Interchangably  In  the  literature  with  minor 
technical  differences  Inserted  by  each  author.  There  are  two  somewhat 
different  types  of  BWERs:  those  of  long  persistence  with  a 5-10  km  dia- 
meter and  those  of  short  duration  with  a diameter  of  only  1-5  km.  A 
dual-doppler  radar  study  of  the  smaller-scale  phenomena  by  McCarthy  et 
il . revealed* dynamics  similar  to  those  discovered  In  the  larger  phenom- 
ena by  Marwitz  and  Berry,  Chisholm  [1970],  and  Foote  et  a^.  [1975];  that 
Is,  both  types  of  regions  contained  a strong  upward  Intrusion  of  warm 
moist  air. 
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CHAPTER  II 
DATA 

Severe  Storm  Date 

unless  otherwise  indicated,  actual  reports  of  severe  storm  occur- 
rences have  been  extracted  from  Storm  Data,  a monthly  publication  of 
the  EnvironmenUl  Data  Service,  National  Oceanic  and  Atmospheric  Ad- 
ministration (NOAA).  This  report  provides  such  information  as  the  lo- 
cation and  time  of  occurrence  of  the  severe-storm  event,  tornado  path 
width  and  length,  maximum  hail  size,  maximum  wind  speed,  and  a brief 
description  of  resultant  injuries  and  property  damage. 

Radar  Characteristics 

The  radar  data  used  in  this  stu(i^'  were  provided  by  the  National 
Severe  Storms  Laboratory  (NSSL),  NOAA.  These  data  were  collected  by  a 
modified  WSR-57  radar  system  developed  and  maintained  by  NSSL  at  Norman, 
Oklahoma.  Table  1 summarizes  the  pertinent  characteristics  of  the  NSSL 
WSR-57  system.  A detailed  description  of  the  design  and  function  of  the 
system  used  for  the  processing  and  recording  of  radar  signals  is  given 
by  Wilk  et  [1967]  and  wm  and  Gra^  [1970].  A block  diagram  de- 
scribing the  acquisition,  processing,  and  recording  of  signals  received 
by  the  NSSL  WSR-57  radar  is  shown  in  Figure  1.  As  shown  in  this  dia- 
gram, the  target  echo  is  detected  and  the  returned  data  processed  gen- 
erally in  the  same  manner  as  a conventional  radar  system.  However, 
after  amplification  and  integration,  the  signal  is  quantized  into. 64. 
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Fig.  1.  Block  diagram  of  the  NSSL  WSR-57  system 
rafter  Wllk  and  Gray.  1970]. 

discrete  Intensity  levels,  the  digital  conversion  process.  The  Inten- 
sity thresholds  for  each  level  are  somewhat  arbitrary  and  vary  slightly 
from  d^y  to  dAy.  In  this  form,  the  data  are  stored  on  magnetic  tapes. 
The  digital  data  are  collected  In  200  Intervals  of  range  called  gates, 
each  approximately  1 km  In  length.  The  precise  length  of  the  gate  Is 
recorded  with  each  set  of  data.  Two  degree  azimuth  Intervals  are  used 
for  data  collection  and  successive  antenna  elevation  angles.  In  steps 
of  1 or  2 deg,  up  to  20  deg  are  used.  The  time  given  for  a particular 
display  Is  the  starting  time  for  the  acquisition  of  the  data  required  to 
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TABLE  1.  Characteristics  of  the  NSSL  WSR-57 
Radar  [after  Wilk  et  al_. , 1967] 


Peak  transmitted  power 

450  kw 

Antenna  gain 

7.079  X 10^ 

Pulse  length 

1.20  X 10^  m 

Beam  width 

2.0  deg 

Wavelength 

10.4  cm 

Pulse  repetition  frequency 

164  sec'^ 

Minimum  detectable  signal 

-no  dbm 

compute  the  display*  The  acquisition  time,  which  Is  a function  of  the 
antenna  rotation  rate,  is  approximately  4 min. 


The  Theoretical  Basis 


The  form  of  the  equation  used  In  this  study  follows  the  derivation 
of  Probert-Jones  [1962]  who  assumed  a more  reliable  beam-shape  In  order 
to  reduce  the  error  Inherent  In  a more  conventional  derivation  [Batten , 


1959]. 

If  the  volume  illuminated  by  the  radar  beam  is  filled  by  the  tar- 
get and  there  Is  no  attenuation  of  the  radar  energy  due  to  preclpiu- 
tlon,  the  average  back-scattered  power,  received  from  a meteorolog 
leal  target  at  range,  r.  Is  given  by 


r ■ Z . 
r ^ e* 


(1) 


where  Is  the  dielectric  constant  used  in  scattering  theory  and  is 
the  equivalent  radar  reflectivity  factor.  The  radar  constant,  C,  is 
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determined  by  the  characteristics  of  the  given  radar  system,  viz., 

w^P.G^he^ 

C . —t  . (2) 

1024A‘^ln2 

« 

where  Is  the  peak  transmitted  power,  6 the  antenna  gain,  h the 
pulse  length,  e the  beam  width,  and  x the  wavelength  of  the  radar. 
From  the  values  In  Table  1,  the  radar  constant  for  the  NSSL  WSR-57  has 
a value  of  1.332  x 10^^  w m”^  or  1.332  x 10"^®  w km^  mm 

For  the  purposes  of  this  study,  the  digital  radar  values  have  been 
converted  to  equivalent  radar  reflectivity  values,  First,  each 

digital  Integer  Is  converted  to  a received  power  value,  P^  (In  units  of 
negative  dbm).  This  Information  Is  provided  on  the  magnetic  tape  via 
the  header  record  generated  prior  to  each  scan.  Next,  (1)  Is  solved  In 


terms  of  Z^,  viz.. 


V' 


or 


(3) 


The  logarithm  of  (3)  is 


log  Zg  ■ log  ♦ 2 log  r - log  CK^. 


(4) 


Using  a value  of  0.9  for  FBattan.  1973]  and  the  value  of  C obtained 
earlier,  we  have 


log  • -9.9. 


(5) 


Vor  a discussion  of  radar  reflectivity  factor,  Z,  equivalent  ra- 
dar reflectivity  factor,  Z^,  and  dbm,  see  Battan  [1973]. 
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Substitution  of  this  value  In  (4)  yields^ 


log  Zg  ■ log  Fj.  + 2 log  r + 9.9.  (6) 

However,  since  the  received  power  values  are  calibrated  In  terms  of 
dbm  Instead  of  watts,  (6)  Is  written 

log  Zg  ■ 0.1  + 2 log  r + 6.9,  (7) 

where 


log  Pj.  (watts)  ■ 0.1  P (dbm)  - 3.  (8) 

Finally,  a digital  value  then  may  be  converted  directly  to  a Z^  value 
by 

Ze(D^)  « 10^®*^  **1  (9) 

where  Is  a digital  value  (0  thru  63),  Is  the  corresponding  aver- 
age received  power  In  dbm  as  provided  In  the  header  record,  r Is  the 

range  to  the  power  return  location  In  km,  and  Z^  Is  In  units  of 

6 -3 
inn  m . 

It  should  be  noted  that  threshold  calibration  values,  as  opposed 
to  logarithmic  mean  values  between  thresholds,  have  been  used  In  the 
digital  conversions  In  this  study.  The  actual  reflectivity  factor 
values  will  be  somewhere  between  the  calculated  values  and  the  next 
higher  threshold  setting. 

^The  radar  actually  obtains  log  rather  than  log  F^.  An  equation 
has  been  obtained  by  Wllk  and  Kessler  [1970]  which  corrects  (6) 
slightly.  The  correction.  In  general,  amounts  to  less  than  1 <fi>. 
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Data  Processing 

Beginning  with  the  1973  storm  season,  NSSL  radically  changed  Its 
data  collection  and  recording  procedures.  Range  Intervals  were  con- 
verted from  nautical  miles  to  kilometers,  the  number  of  digital  Inten- 
sity levels  was  Increased  from  9 to  64,  and  the  magnetic-tape  record- 
ing format  was  modified.  Of  these  changes,  the  latter  caused  by  far 
the  greatest  difficulty.  As  a result  of  the  format  modification,  all 
the  data  no  longer  could  be  extracted  from  the  tapes  by  using  standard 
FORTRAN  "READ"  statements;  In  fact,  none  of  It  could. 

Appendix  A contains  the  record  formats  for  the  1974  and  1975  7- 
track  data  tapes  that  have  been  used  In  this  sf'<dy.  As  Indicated  In 
these  formats,  the  first  record  of  each  PPI  sector  scan,  called  the 
header  record.  Is  recorded  by  using  standard  8-b1t  bytes,  readable  In 
FORTRAN.  The  follow-on  data  records,  containing  two  radlals  of  Infor- 
mation each,  are  recorded  In  6-b1t  bytes  as  only  six  binary  bits  are 
needed  to  record  the  digital  Intensity  Integers  0 thru  63.  These  are 
readable  only  through  bit  manipulation  via  FORTRAN  or  through  an  as- 
sembler language  subroutine.  Further,  the  header  record  Is  readable 
via  FORTRAN  bnly  when  using  the  7-track  TRTCH*T  option,  and  attempts 
to  read  the  subsequent  6-b1t  data  records  with  the  TRTCH«T  option  in 
effect  will  produce  an  abnormal  termination  of  the  program.  Hence,  none' 
of  the  data  can  be  extracted  by  using  standard  FORTRAN  "READ"  statements. 

In  order  to  extract  the  data  from  the  tapes,  a lengthy  and  costly 
program  was  written.  Though  an  assembler  language  subroutine  would  ex- 
ecute faster,  and  hence  be  cheaper  to  run,  a FORTRAN  bit-manipulation 
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program  was  chosen  so  that  an  analyst  without  an  assembler  language 
background  would  be  able  to  understand  and  use  the  program.  The  entire 
fully-commented  program  used  to  read  the  tapes  and  develop  the  presen- 
tations Is  listed  In  Appendix  B.  Statement  numbers  11  thru  83  are  used 
to  extract  data  from  the  tapes. 

Presentations 

It  Is  difficult  to  build  a mental  picture  of  the  three-dimensional 
structure  of  a thunderstorm  from  reflectivity  values  collected  In  a 
spherical  array.  If  the  display  Is  made  at  a constant  antenna  tilt, 
great  distortion  can  occur  In  the  display.  For  example.  If  on  the  4-deg 
antenna-tilt  display,  the  near  edge  of  a particular  Isollne  Is  at  a 
range  of  60  km  and  the  far  edge  of  the  same  Isollne  Is  at  80  km,  the 
difference  In  height  of  the  two  edges  would  be  4577  ft.  The  difference 
becomes  greater  with  an  Increase  In  antenna  elevation  angle  and  range 
from  the  radar.  The  picture  Is  further  complicated  by  the  curvature  of 
the  earth  and  the  refraction  of  the  radar  beam.  It  Is  desirable,  there- 
fore, to  present  reflectivity  values  In  a rectangular  display  (x,  y,  h), 
where  x and  y are  the  rectangular  coordinates  on  a flat  earth  and  h Is  a 
constant  height  above  ground  level. 

Constant  Altitude  Reflectivity  Maps 

Greene  [1971]  developed  a technique  to  synthesize  constant-altitude 
reflectivity  (Z)  maps  or  CAZMs.  Vogel  [1973]  modified  Greene's  techni- 
que and  presented  an  explanation  of  the  procedure  In  excellent  detail. 
The  present  stud|y  employs  this  basic  technique  with  minor  programming 
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CAZM  for  15  kft,  1810  CST,  23  H»y  1974. 
Isopleths  of  reflectivity  In  dbz. 
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KILOMETERS  WEST  OF  NSSL 

Fig.  12.  CAZM  for  50  kft,  1810  CST,  23  May  1974 
Isopleths  of  reflectivity  in  dbz. 
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level  Z as  1 nm®  m"^. 

Vertical  Summation  of  Reflectivities 

Although  CAZMs  describe  the  detailed  structure  of  a storm  system, 
a more  rapid,  less  detailed  method  of  analysis  Is  desirable  for  opera- 
tional applications.  Greene  [1971]  developed  a technique  which  Inte- 
grates the  reflectivity  factor  (In  terms  of  liquid  water)  through  the 
vertical  extent  of  the  storm,  thus  reducing  the  analysis  to  a single 
two-dimensional  display.  Greene  achieved  an  Integration  In  depth  by 
averaging  adjacent  5000-ft  CAZM  reflectivities,  converting  the  average 
reflectivity  to  liquid  water,  and  then  summing  the  ten  resultant  liquid 
water  values  over  a given  x,y  coordinate.  The  resultant  display, 
called  VIL,  was  useful  In  Identifying  centers  of  maximum  reflectivity 
Intensity  and  their  relative  magnitudes. 

Canipe  [1973],  finding  that  VIL  hid  many  Important  features  of  a 
storm,  separated  the  total  VIL  Into  three  distinct  layers  called  Par- 
tial Integrated  Liquid  Water  (PVIL).  The  layers  correspond  roughly  to 
the  source  region  (surface  to  15  kft),  the  active  region  (15  left  to 
35  kft),  and  the  Ice  region  (35  kft  to  50  kft).  This  analytical  tool 
provides  Information  on  where  development  Is  taking  place  In  the  verti- 
cal direction  and,  via  an  x-y  plot  of  the  PVIL  centers,  the  dominant 
tilt  of  the  cell. 

After  evaluation,  another  method  of  presenting  reflectivities  was 
found  to  be  more  economical  In  computer  time  than  Can1pe*s  while  retain- 
ing all  the  desired  storm-structure  features.  This  technique  eliminates 
the  depth-averaging  of  reflectivities  for  Integration  and  the  constant 
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algorithm  conversion  from  reflectivity  factor  to  liquid  water.  Instead, 
the  CAZM  reflectivities  above  a given  x.y  coordinate  are  summed  through 
the  lower,  mid,  and  upper  PVIL  layers  and  the  resultant,  summed-total 
reflectivity  factor  converted  to  dbz  for  map  presentation.  Therefore, 
lower,  mid,  and  upper  Partial  Vertical  Sunined  Reflectivity  (PVSZ)  maps 
replace  the  lower,  mid,  and  upper  PVIL  maps,  respectively.  Figures  13 
thru  15  are  PVSZ  maps  for  1810  CST,  23  4tay  1974. 

In  place  of  the  VIL  presentation,  the  ten  CAZM  reflectivity  values 
(in  dbz)  over  a given  x.y  point  are  summed  in  depth.  The  product, 
called  Total  Vertical  Summed  Reflectivity  (TVSZ),  allows  the  investiga- 
tor to  evaluate  Greene ' s explosive-development  theory  while  reducing 
computer  execution  time  (and  cost).  An  example  of  a TVSZ  map  for  1810 
CST,  23  May  1974,  is  shown  in  Figure  16. 
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CHAPTER  III 


CASE  STUDIES 


The  next  three  sections  contain  radar  case  studies  for  selected 
storms  In  central  Oklahoma  during  1974  and  1975.  In  each  case  study, 
the  environmental  conditions  and  severe-storm  events  are  described 
briefly,  folloMed  by  a storm  structure  analysis  based  on  CAZNs,  and 
PVSZ  and  TVSZ  presentations.  The  significant  features  to  be  brought 
out  are  the  existence  and  location  of  BWERs,  the  tilt  of  the  storm  core 
as  revealed  by  PVSZ  analyses,  the  upper -level  mass  changes  as  provided 
by  the  Upper  PVSZ,  and  the  rapid  Increases  In  TVSZ. 

Case  Study  of  the  Storms  on  23  Nay  1974 
Air  Mass  and  Wind  Structure 

The  synoptic  conditions  In  Oklahoma  for  the  evening  of  23  May  1974, 
that  were  favorable  for  the  formation  of  severe  storms,  are  presented 
In  Figure  17.  A slow-moving  cold  front  had  stretched  across  northern 
Oklahoma  for  several  days  with  a series  of  waves  moving  rapidly  north- 
eastward along  the  front.  On  this  evening,  a wave  located  In  the  Texas 
Panhandle  was  accelerating  the  flow  of  warm,  moist  Gulf  air  northward, 
and  thunderstorms  developed  along  the  front  from  southeast  Missouri  to 
eastern  New  Mexico  as  shown  In  Figure  18.  A mid-level  southwesterly 
jet  provided  the  dry  air  aloft  that  produced  the  atmospheric  Instabil- 
ity needed  for  the  development  of  severe  activity. ' The  winds  aloft  for 
Oklahoma  City  at  1800  CST,  23  Nay  1974,  presented  In  Table  2,  show  this 
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TABLE  2.  Wind  Data  for  Oklahoma  City  at  1800  CST.  23  May  1974. 
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mid-level  Jet  at  9000  ft  with  a maximum  speed  of  14  m sec~\ 

Severe -Stone  Events 

Thunderstorms  producing  various  forms  of  severe  weather  developed 
throughout  the  day  along  the  frontal  zone.  Golf-ball  sin  hall  was  re- 
ported In  early  afternoon  In  north-central  Oklahoma.  As  the  system 
moved  southward,  a funnel  cloud  was  reported  at  Kingfisher,  90  km  north- 
west of  NSSL,  at  1640  CST.  A satellite  cell  southwest  of  the  parent 
near  KlngflsUr,  continued  to  develop  and  produced  a tornado  which 
touched  down  briefly  at  1833  CST  In  the  outskirts  of  Yukon,  40  km  north- 
west of  NSSL.  Later,  a cell  to  the  northwest  of  Oklahoma  City  developed 
to  strong  intensity  and  a funnel  cloud  was  reported  over  Tinker  AFB  at 
1930  CST. 

f 

History  of  the  Yukon  Storm 

During  the  afternoon,  the  squall  line  moved  slowly  to  the  southeast 
with  Individual  cells  moving  to  the  east  at  IS  m sec*\  At  1720  CST,  an 
Intense  cell  (Cell  D)  developed  on  the  southwest  end  of  an  Intense 
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portion  of  the  line.  Within  10  min,  another  less  Intense  cell 
(Cell  A)  developed  30  km  to  the  rear  of  Cell  0 and  the  main  line  of 
activity*  At  1740  CST,  Cells  A and  D began  to  sej^arate  from  the  line 
by  moving  to  the  right  (southeast)*  and  strong  satellites  developed  on 
the  southern  extreme  of  each  parent  cell.  Figure  19*  the  TVSZ  analysis 
sunmary*  clearly  shows  the  movement  of  the  centers  of  maximum  reflec- 
tivity which  correspond  to  each  cell  In  the  system.  The  values  of- 
TVSZ*  displayed  at  5-m1n  Intervals*  reveal  the  relative  Intensity  of 
each  cell. 

As  the  cells  moved  to  the  southeast*  the  tilt  of  Cell  0 changed 
from  easterly  to  southeasterly.  Figure  20*  the  PVSZ  analysis  summary* 
shows  the  tilt  of  Cell  0 as  well  as  that  of  the  other  major  cells.  The 
following  system  Is  used  to  depict  PVSZ  centers  In  a time  series.  The 
symbol*  0*  represents  the  location  of  thd  lower  PVSZ  center*  and  each 
succeeding  dot  is  the  middle  and  upper  PVSZ  centers*  respectively.  An 
excellent  description  of  this  display  method  Is  given  by  Canipe  [1973]. 

At  1800  CST*  a BWER  formed  east  of  satellite  Cell  C.  By  1810  CST* 
the  lower  CAZMs  (Figures  3-5)  and  the  Lower  PVSZ  (Figure  13)  clearly 
showed  a pronounced  hook  echo  Just  south  of  the  BWER;  however*  no  tor- 
nado was  repbrted.  The  hook  existed  through  1820  CST  but  dissipated  by 
1825  CST  as  shown  In  Figure  21.  The  BWER*  however*  remained  very  evi- 
dent throughout  the  time  period  with  the  minimum  reflectivity  steadily 
decreasing  from  33  dbz  to  6 dbz.  Additionally*  the  tilt  of  Cell  0 
changed  further  from  southeasterly  to  southerly*  bringing  It  over  the 
BWER  at  1820  CST.  By  1825  CST*  the  gradient  of  reflectivity  over  a 4- 
km  distance  surrounding  the  BWER  had  reached  a maximum  of  49  dbz 


Fig.  19.'  TVSZ  centers  fbr  the  Yukon  stem  on  23  1974 


Fig.  21.  Lower  PVSZ  map  for  1825  GST.  23  Nay  1974. 

Isopleths  of  reflectivity  In  dbz. 

(55  dbz  near  Cell  C to  6 dbz  In  the  BWER). 

For  unknown  reasons*  the  collection  of  multi -tilt  data  was  termi 
nated  at  1825  GST*  8 min  before  the  tornado  was  reported.  Collection 
was  not  reswMed  until  1845  GST*  and  then  only  below  the  4-deg  eleva- 
tion angle.  By  this  time  the  upper  levels  of  the  storm  had  moved  In- 
side the  maximum  collectable  elevation -angle  of  the  radar  so  that  only 
lower  level  values  of  PVSZ  are  available  at  1845  GST.  These  values* 


hoMver*  do  show  that  a significant  decrease  In  PVSZ  occurred  In  the 
lower  levels  between  1825  GST  (before  the  tornado)  and  1845  GST  (after 


the  tornado).  The  development  of  this  system  is  typical  of  storms  in 
Oklahoma  and  folloMS  an  evolution  almost  identical  to  that  of  the  Ninco 
storm  of  26  April  1969  which  is  described  by  Canipe  [1973]. 

Evaluation 

The  evaluation  section  of  each  case  study  compares  the  significant 

features  of  the  storm  system  with  those  that  would  be  expected  from 

using  Canige's  technique  for  identifying  tornadoes.  For  this  case 

study,  the  BWER  associated  with  the  tornado  was  present  for  at  least 

30  min  before  the  tornado  but  had  dissipated  8 min  after  the  tornado 

was  reported.  Hence,  the  appearance  of  a BWER  is  not  sufficient  evi> 

Q 

dence  to  say  that  a tornado  is  occurring. 

Next,  the  tilt  of  Cell  D moved  over  the  BWER  at  1820  CST.  13  min 
before  the  tornado  was  reported.  An  evaluation  of  the  duration  of  this 
conjuncture  is  impossible,  for  by  1825  CST.  the  upper  reaches  of  both 
Cells  C and  0 were  so  close  to  NSSL  as  to  make  measurements  of  upper- 
level  PVSZ  locations  (and  upper-level  mass  changes)  impossible.  For 
the  same  reason,  the  TVSZ  values  after  1820  CST  are  of  1 value,  as 
only  the  lower  portions  of  these  cells  were  sampled  by  the  radar.  By 
combining  the  above  facts,  we  see  that  Canipe *s  identifying  character- 
istics' are  present  13  min  before  the  reported  touch-down  of  the  tornado. 

If  a BWER  is  defined  simply  as  an  area  of  reduced  reflectivity 
where  the  surrounding  eight  values  of  reflectivity  are  greater  than  the 

I 

central  value  (as  suggested  by  Canine),  other  BWERs  can  be  found  in  this 
storm  system.  At  IBl^CST.  a BWER.  shown  in  Figure  22.  formed  northeast 
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KILOMETERS  NEST  OF  NSSL 


Fig.  22.  Lower  PVSZ  map  for  1810  CST.  23  May  1974 
for  Cell  I.  Isopleths  of  reflectivity  In  dbz. 

of  Cell  H and  southeast  of  Cells  I and  E (see  Figure  19}  and  a hook- 
shaped appendage  formed  south  of  the  BWER.  Further,  the  tilt  of  Cell  E 
at  1810  CST,  though  not  over  the  BWER,  was  to  the  right  of  the  cell 's 
movement  and  to  the  right  of  the  tilt  of  other  cells  In  the  Immediate 
area  (see  Figure  20);  and  a significant  decrease  In  Upper  PVSZ  took 
place  between  1800  CST  and  1810  CST  (28  dbz  to  15  dbz).  These  charac- 
teristics are  Identical  to  those  required  by  Canipe  [1973]  and  quite 
similar  to  those  of  the  Yukon  storm;  however,  there  are  some  noticeable 
differences.  The  core  of  Cell  E,  though  tilted  to  the  right,  was  never 


closer  than  18  km  to  the  BWER  and  the  4-lcm  reflectivity  gradient  was  a 
mere  16  dbz  (48  dbz  to  32  dbz). 

Another  BHER  was  visible  at  1810  CST,  Just  west  of  Cell  H,  and  Is 
shown  In  Figure  23.  Again  the  gradient  Is  small,  only  12  dbz  (45  dbz 
to  33  dbz).  An  out-of-positlon  hook-shaped  appendage  also  Is  visible 
yet  no  tornado  was  reported. 


KILOKTERS  WEST  OF  NSSL 


Lower  PVSZ  map  for  1810  CST,  23  Nay  1974 
Isopleths  of  reflectivity  In  dbz. 

Grid  Interval:  2 km  x 2 km. 


Case  Stucty  of  the  Storms  on  8 June  1974 


Air  Mass  and  Wind  Structure 


The  synoptic  conditions  In  Oklahoma  for  the  late  afternoon  of 
8 June  1974  are  presented  In  Figure  24.  A cold  front,  which  had  re- 
mained stationary  for  several  days,  was  suddenly  pushed  southeastward 

by  flow  around  a large  cut-off  upper  low  over  the  Dakotas.  Thunder- 
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storms  broke  out  in  early  afternoon  from  central  Oklahoma  northward 
through  Kansas  and  northern  Missouri  as  shown  in  Figure  25.  A strong 
low-level  Jetstream  over  southeastern  Oklahoma  was  accelerating  moist 
air  northward  at  20  m sec~\  The  winds  aloft  for  Oklahoma  City  at 
1800  CST,  8 June  1974,  shown  In  Table  3,  show  this  Jet  at  5000  ft  with 
a speed  of  18  m sec~^. 


TABLE  3.  Wind  Data  for  Oklahoma  City  at  1800  CST,  8 June  1974 


tW/K  IM/14  IM/14  IM/17  tOO/li  Z»/lt  lM/14  t4(VI4 


tm/H  nvto  HVt*  <40/17  iM/»  tm/a  tw/a 


Seve re-Storm  Events 


Cells  began  developing  In  the  Oklahoma  City  area  shortly  after  noon 
and  grew  to  extreme  Intensity  within  an  hour.  At  1342  CST,  a tornado  . 
touched  down  on  the  edge  of  the  Will  Rogers  Airport  and  traveled  across 


southwestern  Okalhoma  City.  After  lifting  over  the  downtown  area.  It 
touched  down  again  at  1411  CST  In  the  northeast  suburb  of  Spencer. 
Moments  later,  at  1418  CST,  another  tornado  was  reported  only  5 km 
north  of  Spencer  In  the  conmunlty  of  Jones,  Oklahoma. 

As  the  area  of  Initial  development  moved  northeasterly  toward 
Tulsa,  new  cells  formed  rapidly  on  the  southwest  flank  so  that  a line 
had  formed  by  1500  CST.  At  1535  CST,  a tornado  spawned  by  the  south- 
western-most cell,  touched  down  In  southwest  Oklahoma  City,  2.5  km 
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southwest  of  the  first  one,  but  It  remained  on  the  ground  only  briefly. 
At  1555  CST,  a tornado  was  reported  at  Harrah,  34  km  east  of  Oklahoma 
City  and  39  km  northeast  of  NSSL.  (For  location  of  cities,  see  Figure 
28).  Moments  later,  at  1556  CST,  a tornado  touched  down  8 km  southwest 
of  Orwrlght  (95  km  northeast  of  NSSL)  and  stayed  Intermittently  on  the 
ground  for  89  min  while  covering  a distance  of  85  km.  During  this  time, 
over  200  homes  were  destroyed  and  14  people  were  killed.  Tornadoes  con- 
tinued to  fbrm  east  of  Oklahoma  City  and  south  of  Tulsa  throughout  the 
evening,  the  twentieth  and  last  being  reported  at  2224  CST. 

History  of  the  Drumright  Storm 

At  1510  CST,  a group  of  students  from  the  University  of  Oklahoma 
spotted  a large  funnel  cloud  near  Chandler,  Oklahoma,  60  km  northeast 
of  NSSL.  The  students  watched  the  funnel  for  nearly  10  min  as  It  moved 
to  the  northeast.  Interestingly,  the  nearest  cell  to  this  sighting, 
some  12  km  north,  shows  none  of  Canine's  severe  storm  characteristics. 

At  1530  CST,  the  funnel  was  sighted  again  by  the  State  Police, 
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Figure  27,  the  PVSZ  analysis  summary,  shows  the  tilt  of  Cell  B at 
1530  CST  to  be  almost  vertical.  However,  the  tilt  of  the  cell  changed 
several  times  during  the  next  few  minutes.  In  a cyclonically-rotating 
manner,  the  tilt  changed  from  slightly  easterly,  to  northwesterly,  to 
southerly,  and  back  to  easterly  by  1600. CST. 


KltOMETERS  EAST  OF  NSSL 

Fig.  27.  PVSZ  centers  for  the  Drumrlght  storm  on  8 June  1974. 

Cell  B displayed  no  unusual  movement  as  It  progressed  northeast  at 
16  m sec*^  along  with  other  cells  In  the  system.  Although  Figure  28, 
the  TYSZ  analysis  summary,  shows  that  the  TVSZ  value  of  Cell  B decreased 
steadily  In  magnitude,  the  BWER  and  hooK  echo  remained  visible  In  the 
digital  analyses.  However,  the  4-km  reflectivity  gradient  around  the 
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centers  fbr  the  Onnright  stern  on  8 June  1974 


BWER  decreased  from  36  dbz  at  1530  CST  (see  Figure  26)  to  24  dbz  (41 
dbz  to  17  dbz)  at  1550  CST, as  shown  In  Figure  29. 


Fig.  29.  Lower  PVSZ  map  for  1550  CST,  8 June  1974 
Isopleths  of  reflectivity  In  dbz. 


At  1556  CST,  the  tornado  touched  down  8 km  southwest  of  Drumright. 
The  only  observable  changes  In  the  stone  between  1550  CST  (before  the 
tornado  as  shown  In  Figure  29)  and  1600  CST  (after  the  tornado  as  shown 
In  Figure  30)  were  a further  decrease  In  the  4-l(m  reflectivity  gradient 
to  22  dbz  and  the  development  of  a cell,  called  Cell  Q,  on  the  southern 
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KIUNCTERS  EAST  OF  NSSL 

Fig.  30.  L(Mer  PVSZ  map  for  1600  CST,  8 June  1974. 

Isopleths  of  reflectivity  In  dbz. 

flank  of  Cell  B.  Neither  the  tilt  of  Cell  B nor  that  of  Cell  Q was  over 
the  BMER,  the  nearest  being  Cell  B,  tilted  6 km  to  the  north,  nor  were 
the  tilts  of  the  cells  the  same,  as  Is  shown  In  Figure  27.  Cell  B being 
easterly.  Cell  Q being  southwesterly. 

By  1610  CST.  the  BUER  and  hook  had  dissipated;  however,  another 


BHER  and  hook  had  formed  on  the  southwestern  edge  of  Cell  Q as  shown 
In  Figure  31.  The  tornado  remained  on  the  ground  and  struck  the  town 
of  Olive  at  1620  CST.  At  this  time,  a power  failure  at  NSSL  caused  a 
disruption  In  data  collection  which  made  a further  analysis  of  this 
storm  system  Impossible. 


Lower  PVSZ  map  for  1610  CST,  8 June  1974 
Isopleths  of  reflectivity  In  dbz. 


Evaluation 


For  this  case  stuty,  the  BWER  associated  with  the  tornado  was  pre- 
sent for  at  least  25  min  before  the  tornado.  Although  the  BWER  had 
dissipated  by  1610  CST,  another  one  had  formed  so  that  the  characteris- 
tic was  still  present.  The  tilt  of  Cell  B was  pever  over  the  BWER  but 
was  as  close  as  6 km  at  1600  CST,  shortly  after  the  tornado  touched  do«m 


Most  cells  In  the  central  Oklahoma  area  on  this  afternoon  maintained  a 
rotational  wobbling  tilt  like  that  of  Cell  B,  rather  than  a relatively 
constant  tilt  as  In  the  Yukon  storm  case.  This  wobble  was  noted  by 
Canipe  [1973]  in  several  of  his  case  studies.  Tomadic  activity  was 
Indicated  whenever  the  tilt  of  the  cell  became  vertical  or  was  in  the 
direction  of  the  BUER.  This  tilt  correlation  Is  excellent  In  this 
case«  for  a large  funnel  cloud  was  sighted  at  1530  CST  when  the  tilt 
was  vertical*  and  a tornado  was  reported  at  1600  CST  when  the  tilt  was 
toward  the  BWER.  By  1610  CST,  however.  Cell  B was  no  longer  visible  In 
the  lower  PVSZ  (see  Figure  31)  and  the  tilt  of  Cell  Q was  away  from  the 
newly-formed  BUER.  Burgess  [1976],  In  his  storm  damage  analysis,  does 
Indicate  that  the  tornado  was  still  on  the  ground  at  1610  CST.  No  ex- 
planation of  this  deviation  from  Can1pe*s  characteristics  is  available. 

Upper-level  mass  changes,  revealed  by  upper  PVSZ  changes  (see  Fig- 
ure 27),  indicate  that  a decrease  of  5 dbz  took  place  after  the  funnel 
was  spotted  at  1530  CST.  Another  decrease  of  5 dbz  occurred  as  the  tor- 
nado touched  down,  just  as  predicted  by  Greene  [1971]  and  Canipe  [1973]. 
The  TVSZ  analysis  sumaary.  Figure  28,  shows  a steady  decrease  In  TVSZ 
values.  This  Is  just  the  opposite  of  that  expected  by  Greene.  Thus, 
excluding  the  TVSZ  values.  Canine's  characteristics  were  present  for 
the  funnel  at  1530  CST  and  the  tornado  at  1600  CST,  but  the  character- 
istics were  not  present  at  1610  CST.  Whether  this  was  a brief  excep- 
tion or  a precursor  of  the  lifting  of  the  tornado,  as  It  did  briefly 
around  1615  CST  [Burgess . 1976],  Is  unknown  due  to  the  power  failure 
at  NSSL. 

Just  as  In  the  case  of  the  Yukon  storm,  other  BWERs  were  present. 
Figure  32  shows  one  that  formed  southeast  of  Cell  E (see  Figure  28).  The 
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Fig.  32.  Lower  PVSZ  map  for  1600  CST,  8 June  1974 
for  Cell  E.  Isopleths  of  reflectivity  In  dbz. 

4-laB  reflectivity  gradient  of  21  dbz  Is  just  1 dbz  less  than  the  gra- 
dient shown  In  Figure  30.  In  fact.  If  the  20  dbz  Isopleth  In  Figure  32 
Is  connected  by  the  dots  rather  the  solid  line,  the  hook  and  BWER  look 
almost  Identical  to  those  In  Figure  30.  Further,  the  tilt  of  Cell  E, 
which  rotates  cyclonically  as  Cell  B,  Is  just  6 km  north  of  the  BWER  at 
1600  CST,  the  time  of  Figure  32. 

The  Harrah  Storm 

The  cell  producing  the  Harrah  tornado  at  1555  CST,  moved  Into  full 
radar  view  only  moments  before  at  1550  CST.  'The  storm  appeared  to  have 
none  of  Canipe's  severe-storm  characteristics.  Figure  33,  the  lower 
PVSZ  analysis  for  the  area  at  1600  CST,  shows  neither  a hook  nor  a BWER 
In  the  conventional  position  relative  to  Cell  R and  the  tornado;  how- 
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Fig.  33.  Lower  PVSZ  map  for  1600  CST,  8 June  1974 
for  Harrah  storm.  Isopleths  of  reflectivity  In  dbz. 


ever*  the  20  dbz  Isopleth  near  the  bottom  of  the  figure  may  be  a por- 
tion of  a BUER  or  hook  echo*  the  remainder  of  which  Is  hidden  by  the 
Oklahoma  City  ground  clutter  which  extends  to  Just  south  of  the  map. 

In  any  case*  the  hook-tomado-cell  relationship  would  have  to  be  quite 
different  from  that  seen  In  the  Yukon  and  Orumright  storms. 
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Case  Study  of  the  Storms  on  13  June  1975 
Air  Mass  and  Wind  Structure 

The  synoptic  conditions  In  Oklahoma  for  the  early  evening  of 
13  June  1975  are  presented  In  Figure  34.  Although  no  frontal  zone  was 
presents  a weak  trough  was  pushing  Into  western  Oklahoma  as  a result 
of  an  upper  level  trough  moving  through  the  area.  Thunderstorms,  which 
broke  out  ahead  of  the  trough  In  early  afternoon,  developed  Into  an 
area  covering  most  of  Oklahoma  by  early  evening  as  shown  In  Figure  35. 
The  thunderstorms  generally  were  grouped  In  pockets  of  Intense  activity 
which  moved  very  slowly  to  the  southeast  at  10  m sec'^.  A moderate  low 
level  Jet  to  the  east  of  NSSL  was  carrying  moist  Gulf  air  northward 
while  dry  air  from  New  Mexico  was  carried  eastward  over  the  entire 
area  by  a strong  Jet  at  mid-levels.  The  winds  aloft  for  Oklahoma  City 
at  1800  CST,  13  June  1975,  shown  In  Table  4,  show  this  Jet  near  10  kft 
with  a speed  of  approximately  15  m sec~\ 


TABLE  4. 

Wind 

Data 

for  Oklahoma  City  at  1800  CST, 

13  June 

1975. 

Hi4#«  (kn) 

tfc 

• 

M U 14.  » 

M 

30 

44 

MM  .1 

1«/M 

toi/n 

t7V14  IM/lt  ZWVU  320/10 

334/10 

324/13 

2M/22 

iMr  mz  tytrsfi  mIaO  ■ 220/10 
mo  PVSZ  «Mr*«o  mIaO  • 310/12 

Moor  2VU  AvtrAii  wtoO  ■ 304/20 

Severe  Storm  Events 

Thunderstorms  developed  In  early  afternoon  with  a cells  growing 
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to  severe  Intensity.  At  1430  CST,  golfball-size  hail  was  reported 
at  the  Stillwater  airport,  108  km  north-northeast  of  NSSL.  Hail,  up  to 
2 in. In  diameter,  was  reported  in  north-central  Oklahoma  near  Medford 
at  1700  CST.  The  southwestern-most  cell  in  this  system  continued  to 
develop  and  a tornado  touched  down  at  1733  CST  in  northwest  Stillwater. 
It  traveled  across  the  city,  lifting  some  10  min  later  in  the 
southeastern  outskirts.  The  tornado  damaged  numerous  homes  and  de- 
stroyed 20  trailers  in  a mobile  home  park.  Only  six  people  were  in- 
jured and  no  one  was  killed. 

As  the  area  of  development  moved  southeastward,  a tornado,  which 
damaged  two  homes,  was  reported  near  Kendrick  (90  km  northeast  of  NSSL) 
at  1845  CST.  At  1712  CST,  a cell,  trailing  the  one  near  Kendrick, 
spawned  a tornado  which  touched  down  only  6 km  southeast  of  the  first. 
This  storm  carried  a hayballer  about  1/4  ml  and  destroyed  a concrete 
block  garage.  Isolated  occurrences  of  hall  were  reported  later  In  the 
evening  near  Tulsa. 

History  of  the  Stillwater  Storm 

The  cell  spawning  the  tornado  at  Stillwater,  called  Cell  A,  moved 
Into  full  radar  view  at  1630  CST.  As  the  cell  moved  southeast  at  10  m 
sec~\  as  shown  In  Figure  36,  Its  TVSZ  value  remained  almost  constant, 
as  did  Its  tilt  to  the  south,  shown  In  Figure  37,  the  PVSZ  analysis 
summary  for  Cell  A.  At  1730  CST,  3 rain  before  the  tornado  touched 
down,  a weak  BHER  and  hook  echo  developed  west  of  Cell  A as  shown  In 
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Fig.  37.  PVSZ  centers  for  the  StIllMater  storm  on  13  June  1975. 

Figure  38.  The  4-lan  reflectivity  gradient  around  the  BWER  was  only  18 
dbz  (19  dbz  In  the  BWER  to  37  dbz).  Figure  39  shows  the  BWER-Cell  A 
relationship  5 min  later  (1735  CST)  while  the  tornado  Is  on  the  ground 
In  Stillwater.  By  1740  CST  (3  min  before  the  tornado  lifted),  the  BWER 
had  dissipated. 

Evaluation 

Although  a BWER  was  present  during  the  tornado*  Its  positional 
relationship  was  unlike  any  other  storm  described  previously  (or  any 
described  by  either  Canipe  [1973]  or  Greene  [1971]).  The  classic 
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eastward-moving  cell,  with  the  BWER  to  the  southeast  and  the  tornado  to 
the  southwest,  was  not  present;  Instead,  a super-cell  of  the  type  de- 
scribed by  Marwitz  [1972],  with  the  tornado  located  very  near  the  cell's 
center  (see  Figure  39),  evolved.  While  the  BWER  was  present,  the  tilt 
of  Cell  A was  not  toward  the  BWER  but  almost  90  deg  to  It,  the  closest 
point  being  13  km  to  the  east  at  1730  CST  (see  Figure  37).  Furthermore, 
the  4-lcm  reflectivity  gradient  was  18  dbz,  only  4 dbz  less  than  that  of 
the  false  BWER  Identifier  near  Cell  E In  the  Drumright  storm  (see  Fig- 
ure 32).  There  were  no  noticeable  changes  In  TVSZ  values  and  the  upper 
PVSZ  values,  reflecting  changes  In  mass,  showed  no  marked  decrease  at 
1725  CST  or  1730  CST,  but  rather  a steady  decline  through  1730  CST  and 
then  a sudden  decrease  from  44  dbz  to  40  dbz  as  the  tornado  dissipated. 
No  other  BWERs  were  present  In  any  of  the  other  cells  In  the  system. 

It  Is  apparent  that  Canipe's  technique  for  the  Identification  of  tor- 
nadoes does  not  work  well  with  this  type  of  storm. 
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CHAPTER  IV 


CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 


The  objective  of  this  Investigation  has  been  to  determine  the 
accuracy  and  operational  usefulness  of  a proposed  method  for  tornado 
Identification  from  the  use  of  analyses  of  digital  radar  data.  Three 
cases  Involving  several  tornadoes  were  evaluated  by  using  the  Identifi- 
cation techniques  described  by  Canipe  [1973].  The  stu(1y  led  to  the 
following  conclusions. 

1)  Although  Table  11  of  Canipe 's  study  shows  his  technique  to  be 
100  per  cent  accurate  In  Identifying  tornadoes,  results  from  this  In* 
vestigatlon  show  the  method  to  be  less  than  100  per  cent  accurate. 

Some  of  the  following  specific  conclusions  concerning  Canipe *$  tech- 
nique are  presented  In  support  of  this  statement. 

2)  The  BUER  Is  always  present  when  a tornado  or  funnel  cloud  Is 
reported.  It  generally  forms  before  the  tornado  touches  ground  and 
dissipates  shortly  before  the  tornado  lifts.  Although  no  BUER  was 
found  In  the.Harrah  storm,  the  BUER  may  have  been  hidden  In  the  ground 
clutter  or  Inside  the  radar  coverage. 

3)  BUERs  are  present  not  only  dLr!ng  tornadoes  but  at  other 
times  also.  Further,  the  4-lcm  reflectivity  gradient  surrounding  the 
BUCK,  which  mey  be  used  as  a measure  of  the  BUER  Intensity,  frequently 
Is  the  same  for  cells  that  produce  tornadoes  as  for  those  that  do  not. 
Therefore,  while  neither  the  occurrence  nor  magnitude  of  the  BUER  Is  a 
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positive  Identifier  of  a tornado,  the  occurrence  can  be  used  as  a high 
probability  Identifier. 

4)  The  tilt  of  the  core  of  reflectivity.  If  not  vertical,  1$  not 
alMays  toward  the  BWER,  as  shown  In  the  Stillwater  stoim.  However,  the 
Yukon  and  Drumright  storms  did  have  a tilt  over  or  toward  the  BWER  and, 
therefore,  the  tilt  can  be  used  as  a high -probability  Identifier. 

5)  The  values  of  TVSZ  revealed  only  one  occurrence  of  explosive 
development  and,  therefore,  a poor  correlation  with  tornado  occurrences. 
This  conclusion  was  Just  as  anticipated  by  Vogel  [1973]. 

6)  In  two  out  of  three  cases  (Stillwater  being  the  exception), 
the  upper  PVSZ  values  decreased  as  the  tornado  touched  down.  Just  as 
predicted  by  Canipe  [1973].  Table  5 Is  a summary  of  the  observed  char- 
acteristics of  each  storm. 


TABLE  5.  Summary  of  Observed  Characteristics 
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I Several  conclusions  also  were  reached  concerning  the  use  of 
Canlpe*s  technique  on  an  operational  basis. 
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1)  PVSZ  analyses  do  provide  a rapid,  simple  tool  for  evaluating 
the  structure  and  relative  Intensity  of  cells  In  a storm  system.  The 
method  of  summing  reflectivities  In  depth  and  then  converting  the  total 
to  dbz  for  display  Is  superior  to  previous  methods  In  that  It  retains 
small-scale  features  such  as  the  BWER  and  hides  few  signatures  present 
In  the  Individual  CAZHs.  The  method  also  takes  less  computer-process- 
ing time  (and  hence  less  money)  than  the  liquid-water  conversion  method 
employed  by  Canipe  [1973]. 

2)  A computer  analysis  of  storm  tilt  could  be  very  difficult  to 
achieve  In  complex  storm  situations.  Middle-  and  upper-level  PVSZ  max- 
ima frequently  exhibit  a pattern  quite  different  from  that  found  In  the 
lower  PVSZ  map.  Two  or  three  lower  PVSZ  maxima  appear  to  merge  Into 
one  centrally-located  mid-level  maximum.  This  makes  the  correlation 
of  a mid-level  PVSZ  maximum  to  Its  low-level  counterpart  a difficult 
task,  particularly  for  a computer  and  Its  programmer. 

3)  Since  all  BWERs  do  not  Indicate  the  presence  of  a tornado,  an 
attempt  was  made  to  ascertain  whether  those  associated  with  tornadoes 
exhibited  some  wilque  characteristic.  Possibilities,  such  as  a limit- 
ing maximum  value  of  the  reflectivity  In  the  BWER,  a limiting  minimum 
value  of  the  reflectivity  of  the  nearest  cell,  and  a limiting  minimum 
value  of  the  reflectivity  difference  over  a 4-km  radius  area  around  the 
BWER,  were  examined.  No  conclusion  can  be  drawn  due  to  the  limited  num- 
ber of  case  studies;  however,  for  informational  purposes,  the  maximum 
value  In  the  BWER  appears  to  be  no  greater  than  22  dbz,  the  maximum  re- 
flectivity value  of  the  nearest  cell  appears  to  be  at  least  46  dbz,  and 
the  4-km  reflectivity  gradient  around  the  BWER  appears  to  be  at  least 
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18  dbz.  Further,  extreme  caution  must  be  exercised  when  using  these 
values,  as  BWERs  not  associated  with  tornadoes  approached,  or  In  one 
case  exceeded,  these  limits. 

Recommendations 

The  following  recommendations  are  made  concerning  the  tornado 
Identification  technique. 

1)  Additional  storm  cases  In  Oklahoma  should  be  studied  to  deter- 
mine statistically  the  accuracy  of  this  Identification  technique  and 
how  Its  accuracy  compares  to  other  presently-used  methods.  Including 
doppler  radar. 

2)  The  additional  storm  case  studies  also  should  be  used  to  deter 
mine  whether  the  limiting  values  placed  on  the  BWER  can  be  refined  sta- 
tistical 1y. 

3)  Storm  cases  from  southern  areas  such  as  Texas  and  the  Gulf 
states  should  be  examined  to  determine  the  applicability  of  this  tech- 
nique to  those  locations. 

The  following  recommendations  are  made  concerning  the  processing 
of  digital  radar  data. 

1)  The  mixed  6-b1t  and  8-b1t  word-length  method  for  storing  digi- 
tal radar  data  on  tape,  presently  employed  by  NSSL,  should  be  abandoned 
In  favor  of  the  previously-used  8-b1t  method  due  to  Incurred  Increased 
cost  of  retrieval  from  the  mixed  word-length  tapes  (approximately  three 
times  that  of  the  older  system). 

2)  Action  should  be  taken  by  NSSL  to  reduce  excessive  noise  ap- 
pearing in  the  tapes  which  would  be  difficult  to  handle  In  an  automatic 
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computerized  Identification  technique. 

3)  Radar  range  markers  frequently  appear  In  the  digital  data. 

They  must  be  eliminated  as  the  markers  Induce  erroneously  high  values 
of  reflectivity  along  their  path. 

4)  A simpler  program  requiring  less  computer  core  space  must  be 
developed  before  a mini-computer  can  be  used  to  provide  real-time  dis- 
plays, either  via  a printer  or  a cathode  ray  tube. 

Concluding  Remarks 

It  has  been  the  purpose  of  this  Investigation  to  determine  whether 
tornadoes  produce  a unique  signature  In  analyses  of  digital  radar  data. 
No  unique  Identifying  signature  was  found,  but  rather  a combination  of 
features  which  Indicate  the  presence  of  a tornado  with  a high  probabil- 
ity. Additional  studies  are  required  to  determine  the  precise  proba- 
bility. The  author  believes  that  work  should  continue  In  this  area  and 
that  techniques  should  be  developed  to  make  this  type  of  Identification 
technique  an  operational  reality* 
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APPENDIX  A 

This  appendix  contains  the  format  used  by  NSSL  for  recording 
digital  radar  data  during  the  years  1974  and  1975.  Both  Appendix  A 
and  Appendix  B are  provided  as  an  aid  to  Investigators  who  may  choose 
to  work  at  some  future  time  with  the  type  of  NSSL  data  used  In  this 
thesis.  It  Is  hoped  that  by  listing  the  actual  program  used  In  this 
Investigation,  other  researchers  will  find  this  documentation  adequate 
as  a starting  point  for  their  work. 


MASTER  MAGNETIC  TAPE 
RADAR  DIGITAL  ENCODER 
RECORD  FORMAT 
1974-1975 


FIRST  RECORD  OF  EACH  PPI  SECTOR 

Characters 

Fortran 

Specification 

Legend 

1-6 

16 

Date  (month-day-year) 

7-12 

16 

Time  (CST) 

13-15 

F3.1 

Tilt 

16-21 

F6.3 

Delay  In  kilometers 

22-27 

F6.3 

Gate  length  In  kilometers 

28-30 

13 

Beginning  azimuth 

31-33 

13 

Ending  azimuth 

34-37 

14 

Normalization  range  In  kilometers 
(zero  entry  Indicates  no  normal- 
ization) 

38 

11 

Hardware  normalization  If  one 

39 

11 

Anomalous  propagation* 

40 

11 

Interference* 

41 

11 

Time  constant 

42-44 

13 

Antenna  rotation  time  In  seconds 

45-47 

F3.1 

Horn  noise  (dB) 

51-306 

64F4.1 

Calibrations  (dB) 

307-309 

13 

Number  of  missing  radlals 

310-438 

NI3 

Azimuths  of  missing  radlals 

Data  records  follow 

* Zero  entry  Indicates  undetermined 

Tapes  are  7 track.  556  BPI.  odd 
parity,  binary 
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DATA  RECORDS  OF  EACH  PPI  SECTOR 

Data  records  on  the  1975  archived  WSR-57  radar  data  tapes  are 
438  character  long  and  contain  two  radial s of  Information  each. 

The  format  Is: 


First  Radial 


1-200 

200  Data  Gates  (0-63) 
(6  binary  bits/gate) 

201-203 

Azimuth 

204-205 

Elevation  Angle 

206 

STC  1»on  0-Off 

207-209 

Julian  Day 

210-215 

Time  (CST)  HHMMSS 

216-217 

Delay  Code 

218 

Gate  Length  Code 

219 

Time  Constant  Code 

Second  Radial 

220-419 

200  Data  Gates 

420-422 

Azimuth 

423-424 

Elevation  Angle 

425 

STC  1«on  0-Off 

426-428 

Julian  Day 

429-434 

Time  (CST)  HHHHSS 

435-436 

Delay  Code 

437 

Gate  Length  Code 

438 

Time  Constant  Code 
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Innedlately  following  the  last  radial  In  the  PPI  sector,  a 
special  flag  radial  will  appear  In  either  the  "left”  or  "right  half" 
of  the  current  record.  This  radial  consists  of  200  gates  of  cali- 
bration Information*  and  19  columns  of  eight's  where  the  housekeeping 
would  normally  appear.  For  example.  If  the  last  data  radial  appeared 
In  column  220-438  (right  half)  then  the  next  record  would  look  like: 

col  1-200  Calibration  64F3.1  ** 

201-219  Eight's 
220-438  Nothing 

If  the  last  data  radial  appeared  In  columns  1-219  (left  half),  the 
record  would  look  like: 

col  1-200  Data 

201-219  Housekeeping 
220-419  Calibration  64F3.1  ** 

420-438  Eight's 

Immediately  following  the  last  record  of  the  last  PPI  sector 
on  the  tape,  there  Is  a record  containing  438  nines,  and  then  an 
ENO-OF-FILE. 

* Calibration  adjusted  according  to  Sirmans.  See  NSSL  Tech  Memo  64. 
Does  not  apply  to  1974  tapes. 

**  Calibration  values  are  recorded  MOO  100;  I.e.,  If  actual  value  Is 
104.5,  recorded  value  Is  4.5. 


APPENDIX  B 


This  appendix  contains  the  program  used  to  read  the  tapes, 
produce  the  digital  analysis,  and  print  the  resultant  maps.  The 
program  Is  a highly  modified  version  of  that  used  by  Phillips  [1975] 
and  Is  based  on  the  algorithms  developed  by  Greene  [1971]. 
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